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Abstract. We revisit the construction of puncture black hole initial data in the 
conformal thin-sandwich decomposition of Einstein's constraint equations. It has been 
shown previously that this approach cannot yield quasiequilibrium wormhole data, 
which connect two asymptotically flat spatial infinities. This argument does not apply 
to trumpet data, which connect the spatial infinity in one universe with the future 
timelike infinity of another. As a numerical demonstration we present results for a 
single boosted trumpet-puncture black holes, constructed in the original version of the 
conformal thin-sandwich formalism. 



1. Introduction 

Numerical simulations of black hole spacetimes have experienced a dramatic 
breakthrough (see [U El E] as well as numerous later publications). Many of these 
simulations now adopt some variation of the Baumgarte-Shapiro-Shibata-Nakamura 
formulation jH E] together with the moving-puncture J2J E] method to handle the black 
hole singularities. 

Typically, moving-puncture simulations adopt initial data that are constructed 
using the puncture method [6J. As we explain in more detail below, the central idea of 
the puncture method is to write the unknown functions as sums of analytically known 
background terms that capture the singularities, and correction terms that are unknown 
but regular. If all goes well, the equations for the correction terms can then be solved 
everywhere, without any need for excision or other means of dealing with the black 
hole singularity. Most applications of this method have adopted Schwarzschild data 
on a slice of constant Schwarzschild time as the background solution. In the Penrose 
diagrams of Fig. [T] below, these slices would be represented by lines connecting spatial 
infinity in one universe with spatial infinity of another universe. The resulting initial 
data represent wormhole data. Most applications to date have applied the puncture 
method in the context of the conformal transverse-traceless (CTT) decomposition of 
Einstein's constraint equations [H El El EE] • 
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The moving-puncture method is used in dynamical simulations, and is based on a 
set of empirically found coordinate conditions, namely the "1+log" slicing condition for 
the lapse [UJ and a "f-freezing" gauge condition for the shift [12J. As demonstrated 
by [13j HU [151 EE El) dynamical simulations of a Schwarzschild spacetime using these 
coordinate conditions settle down to a spatial slice that terminates at a non-zero areal 
radius. An embedding diagram of such a slice, which suggests the name trumpet data, 
is shown in Fig. 2 of [16J (see also [TS] IT!?] for an analysis of these slices). We also 
include such a slice in the lower panel of Fig. [TJ which demonstrates that trumpet slices 
connect spatial infinity in one universe with future timelike infinity of another universe. 
Trumpet-puncture initial data, that adopt the puncture method in the CTT formalism 
and produce black holes in a trumpet geometry, have recently been constructed in 

[201 EI]. 

While most initial data that are currently being used in moving-puncture 
simulations have adopted the CTT formalism, the conformal thin-sandwich (CTS) 
formalism [221 El [23], UU] offers an attractive alternative. Several time derivatives 
can be set to zero in the CTS formalism, suggesting that the resulting data may 
represent quasiequilibrium initial data more faithfully than CTT data. Moreover, the 
CTS formalism yields, as part of the solution, the coordinate system in which these time 
derivatives vanish. 

Given these considerations it is of interest, at least as a matter of principle, to 
construct puncture initial data in the CTS formalism. In [21], Hannam, Evans, Cook 
& Baumgarte asked "Can a combination of the conformal thin-sandwich and puncture 
methods yield binary black hole solutions in quasiequilibrium?" and concluded that this 
is impossible, at least for wormhole data that connect two spatial infinities. Here we 
present numerical examples for boosted black holes to demonstrate that it is possible to 
construct quasiequilibrium trumpet-puncture initial data that connect spatial infinity 
in one universe with future timelike infinity of another, at least in the original version of 
the CTS formalism. The extended version of the CTS formalism, which has been more 
widely used in numerical relativity applications, seems to introduce new complications 
in this approach. 

This paper is organized as follows. In Section [2] we discuss decompositions of 
Einstein's equations and introduce the equations of the CTS decomposition. In Section 
[3] we review the puncture method, go over the arguments of [21] to explain why it cannot 
be adopted within the CTS formalism to construct quasiequilibrium wormhole data, and 
argue that the situation is different for trumpet data. We present numerical results for 
boosted black holes constructed in the original CTS decomposition in Section |4} In 
Section [5] we summarize our results, and briefly discuss this approach in the context of 
the extended CTS decomposition, where new difficulties seem to arise. We also include 
a brief Appendix that reviews the Schwarzschild solution in a maximally sliced trumpet 
geometry. 
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2. Basic equations 

In a 3+1 decomposition of general relativity (see, e.g., [2SJ El EE]), the spacetime M is 
foliated by a family of spatial hypersurfaces S. The ten Einstein equations for the 
spacetime metric g a b are projected into the spatial slices, which results in a set of 
constraint equations and a set of evolution equations. The constraint equations constrain 
the gravitational fields on each spatial slice, while the evolution equations govern their 
evolution from one slice to the next. Constructing initial data for Einstein's equations 
entails finding solutions to the constraint equations. 
It is convenient to write the spacetime metric as 

ds 2 = -a 2 dt 2 + 7 ii (dx i + /3 i dt)(dx j + ftdt), (1) 

where a is the lapse, (3 l is the shift vector, and the spatial metric 7^ is the projection of 
the spacetime metric g a b into the spatial slices S. We also define the extrinsic curvature 

as 

Kij = - ^ (d tlij - Di(3j - Djpi) , (2) 

where the operator Di denotes a covariant derivative with respect to 7^. In vacuum, 
Einstein's equations are then equivalent to the two constraint equations 

R + K 2 - KijK ij = (3) 
(the Hamiltonian constraint) and 

D j (K ij -f j K)=0 (4) 
(the momentum constraint), and the evolution equation 



d t K l3 = a(R l3 - 2K lk K k 3 + KK i3 ) - DiD 



+ (3 k d k K i3 + K lk d 3 (5 k + K k3 d$ k . (5) 

Here K = j^Ky is the trace of the extrinsic curvature, also called the mean curvature, 
and Rij is the Ricci tensor associated with the spatial metric. The above set of equations 
is often referred to as the Arnowitt-Deser-Misner, or ADM equations [25] . 

The constraint equations ^ and Q can be decomposed further with the help of a 
conformal transformation 

lij = (6) 

where ip is a conformal factor and 7ij a conformally related metric. We also split the 
extrinsic curvature into its trace K and trace-free part Ay, 

Kij = Aij + - li3 K, (7) 

and conformally transform Ay according to A^ = ip~ 2 Aij. The Hamiltonian constraint 
Q then becomes 

D 2 ip - \^R - ^ 5 K 2 + \.- v .l, r l" = 0, (8) 
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while the momentum constraint (pi) is 

DjA ij - -^fWjK = 0. (9) 
3 

Here Di, D 2 = D l Di and R are the covariant derivative, the covariant Laplace operator 
and the Ricci scalar associated with 7^. Different decompositions of the constraint 
equations proceed with different further decompositions of A^. 

In the conformal transverse-traceless (or CTT) decomposition, A^ is decomposed 
into a transverse and a longitudinal part. The transverse part can then be chosen freely, 
while the longitudinal part has to satisfy the momentum constraint. In addition to the 
transverse part of A^, the mean curvature K and the conformally related metric 7^ are 
freely specifiable and have to be chosen before the equations can be solved. 

The CTT decomposition as the very attractive feature that, for vacuum, maximal 
slicing (K = 0) and conformal flatness (7^ = rjij, where 77^ is the flat metric in 
any coordinate system), the momentum constraint decouples from the Hamiltonian 
constraint and can be solved analytically for both boosted and spinning black holes. 
These Bowen-York solutions (26] can then be inserted into the Hamiltonian constraint, 
which is the only equation that needs to be solved numerically. 

In the conformal thin- sandwich (or CTS) decomposition, we invoke eq. ^ to write 
Aij as 

A ij = ^ ({Lpf - &) , (10) 

where we have introduced the "densitized lapse" a with the conformal rescaling a = ip 6 a, 
and have defined 

Uij = dcfij. (11) 
The longitudinal operator, or vector gradient L is 

(L(3) ij = &fi + D j (5 i - 2 i lj D k (3 k . (12) 

3 

Inserting (10) into the momentum constraint (Joj) yields 

(A L /3) 4 - {Lf3) ij D j ma = aD^aT 1 ^) + \a^ % & K , (13) 

(A L (3y = D^Uf = D 2 ^ + \ ) iril) i V) + Kfl (14) 



where 



is the vector Laplacian, and R l j the Ricci tensor associated with 7^. 



The Hamiltonian constraint rt8j) and the momentum constraint (13) form the basic 
equations of the CTS decomposition in its original form ([22], see also Box 3.2 in [TU]). 
The freely specifiable variables in this decomposition are the conformally related spatial 
metric 7^, its time derivative u^, the mean curvature K, and the densitized lapse a. 



Choices for these functions can be inserted into eqs. ([8]) and (13), which can then be 
solved for the conformal factor ip and the shift vector /3\ Given these solutions, the 
physical solutions for 7^ and can be constructed by reversing the steps above. For 
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a — 1/2 and Uij = 0, the equations of the original CTS formalism reduce to those of 
the CTT formalism (given suitable choices of the free variables in that formalism); in 
particular, the Bowen-York solutions are solutions to the original CTS system under 
these conditions (see [27]). 



In the extended CTS decomposition (see, e.g., [23]), eqs. (pi) and (13) are 
supplemented with the trace of the evolution equation Combining this trace with 
the Hamiltonian constraint (|8| we obtain an equation 



D 2 (p = <P ( V M,,.!" + + ^Rj - i> 5 d t K + ifFDiK. (15) 

for the product of the lapse a and the conformal factor ■0,0 = aip = aip 7 ■ Instead of the 
conformally rescaled lapse a, the time derivative of the mean curvature dtK (together 
with K, and U^) now acts as freely specifiable variable (see also Box 3.3 in |10j). 

3. The puncture method 

Solving any of the above equations for black holes requires handling the singularities in 
the black hole interior. 

One approach is to excise the black hole interior with the help of suitable boundary 
conditions on the black hole horizon (see, e.g., [281 129] in the context of the CTT 
decomposition, and [301 ED E21 EH] for the CTS decomposition). An attractive feature of 
the excision method is that the boundary conditions can be derived from geometrically 
motivated conditions, and can be used to force the black hole to be in equilibrium. 
However, this approach can lead to complicated numerical grids, and can be tedious 
to implement numerically. Moreover, this approach results in valid data only in the 
black hole exterior, while moving-puncture evolution codes need initial data everywhere 
(compare pU [35] ) . 

The puncture method, on the other hand, is much easier to implement numerically. 
The central idea of the puncture method is to write all functions as sums of a background 
term, describing a Schwarzschild black hole, and a correction term (see [6], see also 
[35] I3T]). The background terms are given analytically and absorb the singularities, 
while, if all goes well, the correction terms remain regular everywhere and can be found 
numerically, without the need of excision, on R 3 . 

The Schwarzschild background terms can be given in a number of different 
coordinate systems. Most applications to date have adopted the data on a slice of 
constant Schwarzschild time (e.g. [5J EH]); these data connect spatial infinity in one 
universe with spatial infinity in another universe; the resulting initial data therefore 
represent wormhole data (see the top panel in Fig. [T]). 

It is also possible to adopt maximally-sliced trumpet data as background data in 
the puncture method [20] [21]. These data appear as the limiting case of the family of 
time- independent, maximal slices of the Schwarzschild spacetime (see, e.g., [39 1 PfO l l4"T]. 



see also Appendix A). Trumpet data connect spatial infinity in one universe with future 



timelike infinity of another universe (see the bottom panel in Fig. [TJ see also Fig. 4 
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Figure 1. Penrose diagrams for a Schwarzschild black hole. The dotted lines are lines 
of constant areal radius R; the Killing vectors are tangent to these lines. The 
thick solid line in the top panel represents a spatial "wormhole" slice connecting the 
spatial infinities in two separate universes. The thick solid line in the bottom panel 
represents a spatial "trumpet" slice. This slice starts at the asymptotically flat spatial 
infinity of one universe, enters the black hole region, and asymptotes to the future 
timclikc infinity of the other universe without ever leaving the black hole region. The 
long-dashed line in the lower panel marks the limiting surface at R — 3M/2. 

in [E]), and approach asymptotically a sphere of proper radius R = 3M/2. When 
displayed in an embedding diagram, these slices resemble the shape of a trumpet (see 
Fig. 2 in [IE]), which motivates the name of the resulting geometry. As we discussed 
in the introduction, these data have the appealing feature that they are close to the 
geometry realized in dynamical simulations in the moving-puncture formalism, and 
avoid a coordinate transition from a wormhole to a trumpet geometry. It turns out 
that these data also have very interesting properties from the perspective of the CTS 
formalism. 

To date, the puncture formalism has been adopted almost exclusively in the context 
of the CTT decomposition, where only the Hamiltonian constraint needs to the solved 
numerically (see [1211131 HI] for exceptions). Given the success of the puncture method in 
the CTT decomposition, it is tempting to apply the same method in the CTS formalism. 
In [21] , the authors therefore asked "Can a combination of the conformal thin-sandwich 
and puncture methods yield binary black hole solutions in quasiequilibrium?" and 
concluded that this is impossible, at least for wormhole data. The argument goes as 
follows. 

Consider a Penrose diagram for a Schwarzschild black hole, as in Fig. [lj The top 
panel in the figure includes a spatial wormhole slice, which connects the spatial infinities 
in two separate universes. Any such slice must, at some point, be tangent to a surface of 
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constant areal radius R, marked as the dotted lines in the diagram. This point, marked 
by a dot in Fig. [TJ defines the slice's minimal surface, or its "throat" . Note that the 
Killing vector £^ , which is timelike in the exterior of the black hole but spacelike inside 
the horizon, is tangent to the surfaces of constant R. For equilibrium data, we require 
that the time vector 

t a = an a + f3 a , (16) 

where n a is the normal on the spatial slice, be aligned with the Killing vector 
this defines the "Killing lapse" ax and the "Killing shift" fi K . At the throat, t a must 
therefore be tangent to the slice, which implies that at this point the Killing lapse ax 
must vanish (see, e.g., equation (28) below). 



Unfortunately, the CTS formalism requires dividing by the lapse in equation (10), 
so that the equations cannot be solved without special treatment of the throat. This 
defeats the purpose of the puncture method. The authors of [21] therefore concluded 
that is is impossible to construct quasiequilibrium, wormhole puncture initial data in 
the CTS decomposition. One way to avoid these problems is to relax the assumption of 
quasiequilibrium, and make different choices for the lapse (see |42l 03] ; see also Section 



4.1 below). 

Another way of avoiding these problems is to relax the assumption that the method 
be based on wormhole data that connect the two spatial infinities of two separate 
universes. We now see that trumpet data offer a very attractive alternative in this 
context. Since these slices do not connect two spatial infinities, but rather a spatial 
infinity in one universe with a future timelike infinity of another universe, these slices 
become tangent to the surfaces of constant R only asymptotically (see also [13]). This 
is demonstrated in the lower panel of Fig. [TJ which includes such a trumpet slice. 
Accordingly, the Killing lapse vanishes only at the end of the slice, where singularities 
are absorbed in the "puncture" . This observation suggests to try to construct trumpet- 
puncture initial data in the CTS formalism. 



4. Numerical Results 



In the original CTS decomposition we need to solve the Hamiltonian constraint ([8| and 



the momentum constraint (13) together with expression (10) for the extrinsic curvature. 



Before proceeding we need to make choices for the freely specifiable variables; we will 
assume conformal flatness 7^ = rjij (which implies R\j = 0), maximal slicing K = 0, 
and, consistent with the notion of quasiequilibrium, we will assume ziij = 0. We will 
postpone making choices for a until we specialize to wormhole and puncture data in 
Sections 4.1 and 4.2 below. With these choices, the Hamiltonian constraint rt8|) reduces 



to 



D 2 tp 



1 



if 



-7/1.. Aij 



(17) 



while the momentum constraint (13) becomes 
(A L P) i -(Lp) i W j lna = 0. 
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The momentum constraint decouples from the Hamiltonian constraint and can be solved 
independently of a solution for the conformal factor if). We decorate the densitized lapse 
6lq with a subscript zero as a reminder that this quantity is held fixed in the original 



CTS formalism. Given a solution for j3 l we can compute A 13 from (10), which now 
becomes 



A 11 



2a 



(19) 



To apply the puncture method, we write the unknown variables as the sum of 
some background term and a correction term (but do not linearize in the correction 
term, so that the equations remain valid even in the non-linear regime). In principle, 
any choice for the background terms is possible, but we will restrict our analysis to 
background choices that satisfy equations (17) through (19) above. Towards that 



end, we pick a spatial slice of the Schwarzschild spacetime and adopt the data on 
this slice as background data. We will assume that this slice is a member of the 
family of spherically symmetric, time-independent maximal slices of Schwarzschild (see 
[39] 20j HI]). The geometry of this slice determines the spatial metric and the extrinsic 
curvature. Adopting isotropic spatial coordinates we can identify the background 
conformal factor ?/>o as well as the traceless, conformally related extrinsic curvature 
Aq . Also associated with this slice is a Killing lapse ax and a Killing shift (5 l K which, 
in a dynamical evolution, would leave the background time-independent. 
We now write the conformal factor if> as 



ip = ip + u. 
and the shift /3 l as 

0* = # + 



Inserting (20) into the Hamiltonian constraint (]17|) yields 
D 2 u = --(ipo + uY 



'ilr 7 A A ij 



(20) 



(21) 



(22) 



where we have used the fact that the background conformal factor must satisfy the 
Hamiltonian constraint. We also have 

1 

2a ' 



A ij 



2«o 2«q 



^(Lbr + Al 



(23) 



while inserting (21) into the momentum constraint (18) results in 
(A i 6) i -(L^' J D J lna = 0, 



(24) 



since the background term (3q satisfies the constraint (18) itself. Equations (22) - (24) 



form the equations of the original CTS decomposition in the puncture approach. Given 
choices for ipo, /3q and ao we can now solve these equations for the corrections u and b l . 

As demonstrated by Laguna [H], solutions to the above equations for boosted 
(or spinning) black holes are geometrically equivalent to the Bowen-York solutions [26] 
obtained in the context of the CTT decomposition. This is evident for the choice 
a = 1/2, for which the CTS equations reduce to the corresponding equations in the 
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CTT decomposition and for which the shift /3 1 becomes identical to the Bowen-York 
vector potential; for other choices of a the shift vector takes a form that results in 
ua = 0. 



4-1. Wormhole data 



For wormhole data we adopt the Schwarzschild solution on a slice of constant 
Schwarzschild time, expressed in isotropic coordinates. The background conformal 
factor is 

M 



^o = l + 



2r 



and we have 



A 13 
^0 



0. 



(25) 



(26) 



For the background quantities to satisfy equation (19) we therefore have to adopt the 
Killing shift as our background shift, 



K 



0. 



The Killing lapse of this slice is 

1 - M/(2r) 
aK ~ l + M/(2r)' 

As we discussed in Section [3j this lapse vanishes at the throat at r 



(27) 



(28) 



2Af, so that 



eq. (23) cannot be evaluated there. Noting that both sides of equation (19) vanish for 



the background data, we recognize that we can make other choices for the lapse and 



still have the background quantities satisfy equations (17) through (19). Following [42J, 
we therefore choose the background lapse according to 

1 + c/(2r) 

a ° = l + M/(2r)' (29) 
where c is a constant. For c = — M we evidently recover the Killing lapse, but we 
will choose c = M instead throughout this paper. Given «o we compute the densitized 
background lapse from a = a /ipo- 



With our choice of the background lapse (29 ), the background data satisfy equations 



(17) through (19), but they do not satisfy the trace-free part of the time-evolution 
equation ^ with d t K^ = 0. These data would therefore result in a non-trivial time 
evolution. It is in this sense that these data do not represent quasiequilibrium data. 
To construct quasiequilibrium data we would have to choose the Killing lapse for the 



background. While this is not possible for wormhole data, we will see in Section |4.2 
that this is possible for trumpet data. 

For the above choice of qiq, Laguna [H] found an analytic solution to the momentum 
constraint (24) that carries a linear momentum P l L = (0, 0, Pl) and that, by construction, 
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results in the same extrinsic curvature as the corresponding Bowen-York solution. 
Adopting the notation of |42] this solution is given by 

xzf(r, M, c) 



X 



^ 5(M + 2r) 6 ^' 

n „ vzf(r, M,c) „ . , 

PL = - y ^W P - ^ 

4z 2 f(r,M,c)+g(r,M,c) 
Pl " 20(M + 2r) 6 L ' 

where the functions f(r, M, c) and g(r, M, c) are given by 

f(r, M, c) = (M + c)M 2 + 12(M + c)Mr 

+ 60(M + c)r 2 + 160r 3 (31) 

and 

g(r, M, c) = 5(5M + c)M 4 + 60(5M + c)M 3 r 
+ 2(749M+ 149c)M 2 r 2 
+ 8(497M + 97c) Mr 3 

+ 120(49M + 9c)r 4 + 4480r 5 (32) 

(see also the red surface in Fig. [4] below). This solution serves as a very valuable test 
for our numerical code. 

In our finite-difference code, we use Cactus and PETSc software to invert the 



operator on the left-hand side of eq. (24) simultaneously for all three components 
of b l . We adopt Cartesian coordinates with N equidistant gridpoints in each spatial 
dimension. At the outer boundaries, imposed at identical locations A out = F out = Z out , 
we enforce the fall-off conditions of the Laguna solution (30). Wormhole data satisfy 
certain symmetry conditions, so that we could restrict our code to one octant only and 
impose symmetry conditions on the coordinate planes. However, for the trumpet data 
in Section 4J2 these conditions to not apply; we therefore will not use these symmetry 
conditions here either. 



Since eq. (24) is linear, solutions cannot be unique. We pick a particular solution 
by fixing the components of h l at the puncture r = 0, where eq. (24) becomes singular 

b* - 5M + V (33) 

°punc - 

(see |42j for a more detailed discussion and motivation for this approach). We choose a 
vertex-centered grid, so that the puncture coincides with a grid point, and we force the 
solution to take the value (33) there. 

As an example, we show some numerical results for the x component of the shift in 
Fig. [2j demonstrating that our numerical solutions converge to the analytical solution. 
In agreement with the findings of |42], the rate of convergence is between first and 
second order; the lack of second-order convergence is caused by the singular nature of 
the puncture, where we force the solution to take fixed values. 
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Figure 2. Convergence of the shift component f3 x to the Laguna solution (30) for a 
boost of P = 0.5M and c = M. Here the outer boundaries were imposed at X out = AM, 
and we show results along the space diagonal for three resolutions. The upper panel 
shows the convergence of the numerical solutions to the analytical result, shown as the 
solid line, and the lower panel shows the differences between the analytical solution 
and the numerical results. The insets show enlargements of the regions marked by the 
boxes. 



Given a solution for b l , we can compute from (23) and insert the result into the 



Hamiltonian constraint (22). For wormhole data, the right-hand side approaches zero 



at the puncture, so that the equation can be solved without special treatment of the 
puncture. The Hamiltonian constraint is non-linear in u and cannot be solved directly. 
We linearize the equation in corrections to u and solve the resulting equation iteratively 
until the L2 norm of the residual has dropped below a certain tolerance (here chosen to 
be 10~ 9 in code units). We show results for a boost of P = 0.5M in Fig. 



4-2. Trumpet data 

For trumpet data we adopt the Schwarzschild solution on a trumpet slice as background 
data. More specifically, we adopt the limiting case of the family of maximal slices that 
we described in Section [3j Maximally sliced trumpet data are discussed in [14J , and an 
analytic solution in isotropic coordinates is given in [18J. We summarize this solution 
Appendix A 



in 



The background conformal factor tp Q and the background extrinsic curvature A l i 



are given by (A.l) and (A. 10). As before we adopt the Killing shift, now given by 
(A. 5), as the background shift Since neither Aq nor /3q vanish everywhere, we 



now have to adopt the Killing lapse (A. 4) as the background lapse for the background 



quantities to satisfy eq. (19). Unlike for wormhole data, this does not pose a problem, 
since this Killing lapse remains positive for all r > (see also Fig. 2 in [18J). Adopting 
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Figure 3. Numerical results for the corrections u, b x and b z , for both wormhole 
and trumpet data in the original CTS formalism, for Pl — 0.5M. Results are shown 
along the space diagonal, with outer boundaries imposed at X out = 8M and N = 129 
gridpoints in each dimension. 



the Killing lapse as the background lapse will lead to quasiequilibrium data, which 



represents a potential advantage over the approach for wormhole data in Section 4.1 



As we discussed in Section [3j this observation motivates the study of trumpet-puncture 
initial data in the CTS formalism. 



We can now solve equations (22) - (24) with trumpet background data. Unlike for 
wormhole data, the right hand side of the Hamiltonian constraint (22) does not remain 
finite at the puncture for a trumpet background. For non-spinning black holes this 
implies that the correction u has to vanish at the puncture itself (see [201 EI])- Instead 
of solving the Hamiltonian constraint everywhere, we therefore fix the solution to zero 
at the puncture, similar to how we force the shift to take certain values at the puncture. 

In Fig. [3] we show results for the corrections to the conformal factor and the shift, 
and compare with the corresponding results for wormhole data. To specify the solution 
for the shift, we again set the shift correction b l to (33) at the puncture - the wormhole 
and trumpet data therefore agree at the puncture. Since the Laguna solution (30) is not 
a solution to equation (24) with trumpet background data, the two solutions are not 
identical away from the puncture. However, as shown in Fig. [3j the difference between 
the two solutions is not dramatic. This is also evident in the surface plot in Fig. |4j 
where we show the x component of b l , for both wormhole and puncture data, in the x-z 
plane. 

For trumpet data, the linear momentum is no longer given by the parameter P l L in 
(33). We instead compute the linear momentum from 

= JL j A ij d 2 Sj, (34) 
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Figure 4. A comparison of the ^-component for wormhole and trumpet data, for a 
boost of Pl — 0.5M. The blue data (on top on the right and left) represent trumpet 
data, while the red data (on top in the front and back) represent wormhole data. 



and evaluate the integral on the outer boundary of our computational domain. Given 



that the momentum constraint (24) is linear, we now find the linear relation 

P { = 1.43P{ 
for our choice of c = M. 



(35) 



5. Discussion 



We reconsider the construction of quasiequilibrium puncture black holes in the context of 
the CTS decomposition. It has been shown previously that this approach cannot be used 
for wormhole data, which connect the two spatial infinities in two separate universes (see 
|24j). For such data, the Killing lapse has to vanish somewhere away from the puncture 
even for Schwarzschild black holes, so that the CTS formalism cannot be applied (this 
problem can be avoided by choosing a different lapse function, which, however, does not 
lead to quasiequilibrium data, see [42, 43J). Here we point out that this argument does 
not apply to trumpet data, which connect the spatial infinity in one universe with the 
future timelike infinity of another. For these data, the Killing lapse vanishes only at the 
puncture itself, which suggests that it may be possible to construct quasiequilibrium 
trumpet-puncture data in the CTS formalism. 

We verify that it is possible to construct such quasiequilibrium trumpet-puncture 
data in the original CTS formalism, and present numerical results for single boosted 
black holes. We also compare these trumpet data with wormhole data that have been 
constructed using the same formalism, but with a non-Killing lapse. 

We also attempted to solve the equations of the extended CTS formalism for both 
wormhole and puncture data. The equations appear similar to those of the original 



CTS formalism, except that we now solve the lapse equation (15) together with the 



other equations. This, however, introduces an important difference; instead of fixing 
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a, as in the original CTS formalism, we now fix dtK and compute a from if) and 



Inserting this dependency explicitly in equations (J8J) and (15) changes the sign of the 
exponents of if) and in some of the terms. As a consequence, the maximum principle 
can no longer be applied to establish uniqueness of solutions (see [15J US]). 

We have solved the equations of the extended CTS formalism, for a wormhole 
background with a non-Killing lapse, to reproduce the results of [12] ■ These solutions 
exist only up to a certain critical value of the momentum (compare Fig. 8 in [12]), and 
may not be unique (compare [17J US]). For trumpet data, however, we have been unable 
to find any converging solutions, even for very small values of the momentum. While 
we cannot rule out that the fault lies with our numerical code, we suspect that these 
solutions do not exist. The equations form a complicated system of non-linear equations, 
especially for a trumpet background, and we will postpone a further investigation of 
whether these equations allow any regular non-trivial solutions to a future study. 

However, even the trumpet data constructed within the framework of the original 
CTS formalism are an interesting alternative to more traditional black hole initial data; 
as trumpet data they are closer to the geometry realized in dynamical moving-puncture 
simulations than wormhole data, and as CTS data they are represented in a coordinate 
system that is well suited for numerical evolution. 
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Appendix A. The Schwarzschild solution in a maximally sliced trumpet 
geometry 

A family of maximal slices of the Schwarzschild spacetime can be given analytically 
in terms of the areal radius R (see [39l HOj HI]). The limiting member of this family, 
which describes a trumpet geometry, can also be transformed analytically into isotropic 
coordinates, even though the solution can only be given parametrically, with R as 
the parameter [18J. In our numerical implementation, we determine all variables as 
a function of r by creating tables and using suitable interpolation. 
In particular, the conformal factor is given by 



if) 



AR 



2R + M+ (AR 2 + AMR + 3M 2 ) 1 / 2 

~8R + 6M + 3(8R 2 + 8MR + 6M 2 ) 1 / 21 1/2 ^ 
(4 + 3v / 2)(2i?-3M) 

where the isotropic radius r is 

'2R + M+ (AR 2 + AMR + 3M 2 ) 1 / 2 



n l/2 

(A.l) 



(A.2) 
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x 



(4 + 3\/2)(2i?-3M) 



8R + 6M + 3{8R 2 + 8MR + 6M 2 ) 1 / 2 
Asymptotically, ipo behaves as 

'3M\ 1/2 



1/V2 



2r 



1 + 



M 
2r 



r ->■ 0, 



r — > 00. 



(A.3) 



The limit surface r — > corresponds to a sphere of areal radius i? = T/> 2 r — > 3M/2. 
The Killing lapse and shift are given by 

r4\ 1/2 



Ok 



2M 27.1/ 
1 - — + 



R 16i? 4 



and 



/3 



3V/3M 2 



K 



R r 



To determine the behavior of ax in the neighborhood of the puncture, we write 

„ 3M 
R=^ + P 



and expand both (A. 2) and (A. 4) to leading order in p. This results in 

■2k 



r ~ kM 



9 M 



where we have abbreviated n = (4 + 3a/2)/4, and 
2V2 p 



OCR 



3 M' 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



(Ai 



Combining eqs. (A.7) and (A. 8) we can eliminate p and find that, to leading order, ot K 
scales with r , 



3V2 fl 



k KkMJ 

Finally, the conformally related, trace-free part of the extrinsic curvature is 
3\/3M 2 



A* 3 



4r 3 



-(f 3 -3nV), 



(A.9) 



(A.10) 



where n % = x l /r is the spatial normal vector pointing away from the puncture at r = 0. 
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